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Abstract 

We derive the anharmonic correction to the hyperon energy in the bound 
state version of the topological soliton model for the hyperons, and show 
that it represents a negative correction of at most 10% to the energy of the 
bound heavy flavour two-meson system in the case of cascade hyperons. The 
main anharmonic correction arises from the mass term in the Lagrangian 
density. For large meson masses the consistency of the model requires that 
the anharmonic correction decreases as the inverse square root of the mass 
of the heavy flavour meson. 







1. Introduction 



The bound state version of the topological soliton model for the hyperons 
originally proposed by Callan and Klebanov [1] has been shown to provide 
a remarkably good description of the static properties of the hyperons in all 
heavy flavour generations [2-4] . The version of this model in which the soli- 
ton is stabilized by explicit vector meson fields [3] has been shown to provide 
a smooth interpolating model, which possesses chiral symmetry in the light 
flavour sectors and heavy quark symmetry [6] in the limit when the heavy 
meson mass grows beyond limit [7,8]. As the model thus has all the symme- 
tries that are believed to be important in the low energy sector, it should be 
viewed as realistic enough to have predictive power - as e.g. in its prediction 
of the existence of nonstrange pentaquark states [9] . 

The original version of the bound state soliton model [1] is based on the 
direct extension of Skyrme's Lagrangian density [10] to 577(3), and an ex- 
pansion to second order of the isodoublet kaon field which is treated as a 
quantum fluctuation. The resulting harmonic approximation leads to a wave 
equation for the kaon field, which has two bound states, the lower one of 
which describes the A — X system (P-state) and the excited one the A(1405) 
(5-state) [1,2,11]. The cascade and Q~ particles are then constructed as mul- 
timeson states, which in the harmonic approximation do not interact. In this 
paper we investigate the validity of this harmonic approximation by carrying 
out the expansion in the meson field to quartic order. Because of the compli- 
cated form of the stabilizing term in the Lagrangian density of the Skyrme 
model [10] we employ the simpler alternate form proposed by Pari [12]. This 
alternative stabilizing term, which is quartic in the derivatives, is equivalent 
to Skyrme's quartic stabilizing term at the level of 577(2), but not in the 
case of SU(3) fields. In the case of the kaons it has been shown that lin- 
ear combinations of the two alternate quartic stabilizing terms lead to good 
predictions for the strange hyperon energies, although if the alternate (Pari) 
term is employed alone bound rjN states would appear, in contradiction with 
empirical findings [13]. As the kaon bound state energy is fairly insensitive 
to the choice of quartic stabilizing term we shall however here be content to 
consider only the simpler alternate quartic stabilizing term, as this should 
suffice for the purpose of estimating the anharmonic correction. 
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The anharmonic term of lowest order takes the form of a meson-meson 
interaction, which we treat in lowest order perturbation theory. The matrix 
elements of this interaction for two-meson hyperon states - i.e. the cascade 
particles and their heavy flavour analogs - represent the perturbative esti- 
mate of the anharmonic correction to the hyperon mass values, the quality 
of which should be comparable to the standard perturbation theory estimate 
of the electron-electron interaction energy in the case of the He-atom. We 
find that in the case of the cascade type hyperons the anharmonic correc- 
tion amounts to a negative correction of at most 10% to the energy of the 
bound two meson system. In the case of the strange cascade hyperons the 
correction is -26 MeV, in the case of the S c -233 MeV and in the case of 
the predicted hyperons -639 MeV when standard values for the parame- 
ters of the model are used. We show that in the case of very large meson 
mass values the consistency of the model requires that the anharmonic cor- 
rection is roughly proportional to the inverse square root of the heavy meson 
mass. The relative smallness of the anharmonic correction indicates that 
the harmonic approximation is reliable and that the satisfactory predictions 
that it gives for the hyperon spectra and magnetic moments are fairly robust. 

This paper is organized into sections. In section 2 we derive the general 
form for the lowest anharmonic term in the Lagrangian density of the meson- 
soliton system. In section 3 we show that the matrix element of this term in 
the case of a state with two mesons in the ground state can be reduced to 
a fairly simple form. In section 4 we derive an approximate estimate for the 
dependence on the heavy flavour meson mass of the anharmonic correction. 
In section 5 we present the numerical results for the anharmonic correction 
to the cascade-type hyperons. Section 6 contains a summarizing discussion. 



2. The anharmonic Hamiltonian density 
2.1 The harmonic Lagrangian 

In its simplest version the bound state model for the hyperons is based on 
a Lagrangian density for an SU (3)-valued field U that has three components: 
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C = C a + C s + Ccsb- (2.1) 
The first of these is the Lagrangian density of the non-linear a-model: 

= -\flTrL,L^ (2.2) 
where the "left current" U 1 is defined as 

L m = u^d^U. (2.3) 

The second term Cs is a stabilizing term for which a variety of forms has 
been considered [10,12,14]. We shall here employ the form 

£s = ^{{TrL.Kf - (TrL M L M ) 2 }, (2.4) 

that was suggested by Pari [12] as a simpler, although SU(2) equivalent, 
alternative to Skyrme's original quartic stabilizing term [10]. The third term 
£csb takes into account the chiral symmetry breaking meson mass term and 
the S77(3) symmetry breaking implied by the different values for the decay 
constants of the mesons in the different flavour generations. For this term 
we employ the form [15] 

C CSB = S^f^Tr { U + U^-2 } 

+ /> '~ /2m2 Tr{ v / 3A 8 (t/ + E/t)} 

_fl_H Tr{ ^ _ ^x 8 )(Ud^d"U + rfdJJ&rf)}, (2.5) 

where m represents the mass of the heavy flavour meson (K, D or B) and / 
the corresponding decay constant. 

This Lagrangian density has to be augmented by the Wess-Zumino action 

Swz = / d'xe^TriL.L^LpL,}, (2.6) 

which contributes to the energy when the field manifold is larger than 577(2), 
although it cannot be reduced to the form of a Lagrangian density [1,2]. Here 
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Nc is the colour number. 



For the SU(3) field U we adopt the form [1] 

U = y/u^U m y/U^, (2.7) 

where U n is the soliton field 

In this expression u is the usual SU (2) hedgehog field: 

u = e ^ r \ (2.9) 

where 6 is the chiral angle. 

The heavy meson field U m has the form 

^ m = e#j( ^ t M q )}• (2-10) 

Here M represents one of the meson isodoublets (K + ,K°) T , (D°,D~) T or 
(5+B°) T [4]. 

The bound state model for the hyperons is obtained by expanding the 
Lagrangian density (2.1) to second order in the meson field M. This har- 
monic approximation leads to a linear differential equation of second order 
for the meson field M, in the potential field of the S'L r (2)-soliton. The bound 
state solutions to this wave equation represent stable hyperon states. 

In the harmonic approximation the interactions between the mesons in 
the hyperons, that are formed as multimeson states as the cascade and £l~ 
particles, are neglected. We here carry out the expansion in M to quartic 
order in order to obtain the lowest anharmonic correction to the Hamilto- 
nian density for the meson field M. This takes the form of an effective 
meson-meson interaction, which we treat as a first order perturbation. The 
employment of the alternate quartic stabilizing term (2.2) is motivated by 
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the fact that it leads to a much simpler form for the anharmonic correc- 
tion than Skyrme's original stabilizing term. We show that the alternate 
term leads to binding energies for the heavy isodoublet pseudoscalar mesons, 
which are very close to those obtained with the Skyrme term, and hence it 
should suffice for the purpose of estimating the magnitude of the anharmonic 
correction. If a realistic description of the interaction between the soliton and 
the complete pseudoscalar nonet is required one has in general to consider a 
linear combination of the two quartic stabilizing terms [13]. 

In the expansion of the Lagrangian density (2.1) in powers of the meson 
field M the lowest (zero) order term is the usual Lagrangian density of the 
SU{2) soliton [16]. Because of the stability of the soliton only terms of even 
order in M appear. To second order one obtains in the case of an unrotated 
soliton field [1,7] the following expressions: 

£» = -La„M^M + -Lm»{(2^cos0 4- 

_£!!^/? f .£ } M, (2.H.) 

+ (0' 2 - ^-)r ■ VM» r ■ VM 

+M\2cos8 y—(6 2 + — — ) + —sin 9 

-2— -J-{e 2 + — )f • L]M}, (2.116) 

4 2 L = -m 2 MtM + ( X 2 - 1)4 2 ). (2.11c) 

Here the factor \ is the ratio f / f n . Finally the Wess-Zumino term (2.6) gives 
rise to the following contribution to the meson Lagrangian in the harmonic 
approximation [2]: 

N 

Cwz = ij^B^D^M - (^M)tM). (2.12) 
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Here is the anomalous baryon current of the soliton [16], the time com- 
ponent of which has the explicit form 

B° = --±—sin 2 6 6\ (2.13) 

2 r n jL r A 

and Dn is the "covariant" derivative [1] 

D, = d, + ^(Md^+^d,M)- (2-14) 

The second order Lagrangian densities (2.11), (2.12) do, upon a mode de- 
composition of the kaon field, lead to the following wave equation for the 
bound meson: 

a(r)V 2 k + b(r)f ■ Vk - c(r)L 2 k 

-[v (r) + v IL (r)I- L]k - m 2 k + 2u\(r)k 

+d(r)u 2 k = 0. (2.15) 

Here u is the energy of the bound meson {uj < m) and k(r) is the meson 
wave function. The effective (iso)spin operator / for the system is defined as 
I = t/2. The radial functions a, ...A in eqn. (2.15) are defined as 

< \ 2 sin 2 9 , n „ „ . 

a(0 = l + ^ — , (2.16a) 

, . . 4 sinO , sinO „ . 

6(r) = -__(— - - «„6>6>'), (2.166) 



1 //i'2 sin 2 ^ 

^2^*2,^2 ,y<2 



C W = -T7^(^--^), (2.16c) 



^) = l + ^' 2 + 2^), (2.16d) 

ft n sin 2 6/2 6' 2 1 r ,m 2 fl 
«o(r) = -20086—^ - — - 2 ^ 



A sin 2 6/2 4 sin 2 Ql2, n , sin 2 9, 
vidr) = 4^^ + — —^(fl 2 + — ), (2.16/) 
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N c sin?Off 

A ( r ) = - 87r2/2 r2 • ( 2 - 16 ^) 

The ground state solution to the meson wave equation (2.15) is a P-state. 
In the case of the kaon we obtain the ground state energy to be 195 MeV 
when fx = 1.23 f n , and the parameter values f n and e are chosen as in ref. 
[16] so that the nucleon and A33 resonance have their empirical mass values 
(A — 64.5 MeV, e=5.45). This value for the energy of the bound K meson 
leads to satisfactory predictions for the masses of the stable strange hyperons 
[4]- 

In the case of the heavier flavour hyperons the alternate quartic stabilizing 
term (2.11b) does not lead to as satisfactory predictions for the energies of the 
charm and bottom hyperons if the decay constant ratios /d//tt an d fslf-n 
are varied within the presently accepted uncertainty range. Unless these 
decay constant ratios are given fairly large values the D and B mesons will 
be predicted to be overbound by a significant amount. We shall here use the 
value /d//tt = 2.4, which gives the value ojd = 1226 MeV, which is about 100 
MeV too small. If the value for fo/f-n is reduced to the value 1.8 at the upper 
end of the empirical uncertainty range the predicted value for ujd drops to 
only 996 MeV, which is unrealistically low. In the case of the B meson we 
obtain cjb=3528 MeV with the value /s// 7r =2.8, which is somewhat above 
the presently accepted uncertainty range. This value for uj b is about 1 GeV 
below the more realistic value that is obtained with the original Skyrme 
model stabilizing term [4] , which reinforces the conclusion that the alternate 
quartic stabilizing term should not be used for phenomenological predictions 
except in combination with Skyrme's stabilizing term [13]. We nevertheless 
assume that it is sufficiently realistic to permit an estimate of the relative 
importance of the anharmonic corrections to the energy. 



2.2 The anharmonic Hamiltonian density 

We now turn to the derivation of the lowest order anharmonic correction 
to the Lagrangian density (2.11), (2.12) for the meson field by expanding the 
full Lagrangian density (2.1) to quartic order in the kaon fields. The quartic 
anharmonic correction that arises from the non-linear a-model term (2.2) 
has the following form: 



7 



4 4 ) = _L_{,9 M M t M^M t M + M^d^MM^d^M 
■-M Mi),,.]! i) 1 ' M - M^d^Md^M^M} 

1 



-{(i.U /.U • MVM + .1/ ( / • f + f • l)MM^M 



12x 2 / 2 

+^^^[1^+ f) • VM - VM f • (f+ r)M] 

+3M f (f+ r)M ■ [M f VM - VM f M]}. (2.17) 

— * 

Here / and r are defined as the left and right currents for the square root of 
the soliton field: 

1= \fu$3yJu~ v , (2.18a) 



r = ^/U n V\/ul (2.186) 
The explicit forms of these are 

r , f sin9 _ sin 2 9/2^ „ .& sinO.^ , n „ , 

/ = «{ r + -rxf+ )f-ff}, (2.19a) 

1 2r r V 2 2r J ' V ; 

, r sin9^ sin 2 9/2^ „ ,0' sinO.^ 
f=-«{^— r —fxf + (--^—)T-rr}. (2.196) 

The quartic anharmonic correction that arises from the stabilizing term 
C s (2.4) can be expressed in the following form: 

„C4i 1 r,„'9 sin 2 9. t . » 9 sin 2 9.^ , , 
4 4) = ^{(^ 2 + — - - — )r m r B A nn } 

+ iW {5 ^^ " (Bwi)2} - (2,20) 

Here we have introduced two tensors A M „ and which are defined as fol- 
lows: 

V = \{M^M{d^d u M + d u M^d^M) 
• i)„MWM <),,.]! + i),,MMM i),M 
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-M^d u MM ] d^M - MD^MM i)„M 

i^MMi^MM - <),M Mi),, MM 

^^ u cos9sin 2 9/2 . .9' 2 cos9sin 2 9/2, A 
-2M f [ + (— -^JvjMMtM 

„„, + r sm0 ,9' sin9._ r sm0 ,0' sin9._ _ 

-6Mt[— r ro + (- - — )f • ff m ]MMt[— r n + (- - — )f • fr n ]M 

+6M t K(r)(f x f) m MM ] «(r)(f x r) n M 

-iM^n(r)(f x r) m V n MNlM 

+iV n M^K(r)(Txf) m MMM 

-iM^ K(r)(f x f) n V m MMM 

+iV m M^K(r)(Txf) n MMM 
+3iM^K(r)(r x f) m M(V n MM - M^V n M) 
+3M t «(r)(f x f) n M(V m MM - M t V m M)}, (2.21) 

and 

£v = -2(9 M Mt9,M + d v M^d^M) 

+ + r rn l n )M 

+2w(r)[M t (r x f) m d v M + M t (rx f) n d^M 

-d^M\r x r) n M - S^M^f x r) m M]. (2.22) 

The indices m and n in the expression for i? M „ indicate the spatial components 
of 5. Above we have used the notation 

, . sin 2 9/2 
«(r) = (2.23) 

Finally the charge symmetry breaking Lagrangian (2.5) gives rise to the 
following quartic anharmonic term: 

Cffsa = I^mMmM 

L(i _ ^{-K^Kd^d^K + -duK^ KK^d^K 

2/ 2 X 2 3 ^ 3 M 



9 



—-K^dnKK^d^K - -d u K^Kd»K^K 
3 3 M 

+ ^KK\l m r m + r m l m )K 

~K^K[{V m K\l m + r m )K - K\l m + r m )V m K)] 

+ (K t V m K - V m K^K)K\l m + r m )K 

+2KH m KK^r m K}. (2.24) 

The quartic anharmonic Hamiltonian density is obtained by introduction 
of conjugate momentum fields as 

n = ^£-, nt = ^, (2.25) 

and carrying out the Legendre transformation 

W = n f M + M+n - £. (2.26) 

The resulting expression for the anharmonic term in 7i is identical to that of 
the quartic lagrangian densities (2.17), (2.20) and (2.24), with exception for 
an overall sign change of those terms that do not contain time derivatives of 
the meson field or its Hermitean conjugate. 

The quartic Hamiltonian density has the form of meson-meson interac- 
tion. In general the interacting mesons can be in different states, although 
we shall here consider the case of both mesons in the ground state. As the 
ground state is a P-state it is described by a wave function of the form 

M = A^lr-r, (2.27) 



where k(r) denotes the radial wave function and A is the rotation operator, 
which transforms the meson isospin operator into an effective spin operator. 
As the soliton field is rotated by the same SU(2) rotation operator, i.e. 
u — > AuA^, the isospin operator r in the Lagrangian densities is converted 
into the spin operator of the meson as well. Note that the normalization 
condition for the meson wave function k(r) is 

roc 

2 / drr 2 k 2 (r)[u)d{r) + A(r)l = 1, (2.28) 
Jo 
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a result that is implied by the form of the harmonic Lagrangian for the meson. 



The anharmonic correction to the energy of a two-meson state, which 
corresponds to a hyperon of the cascade type, is obtained as 



where H 4 is the Hamiltonian density of quartic order in M. The expression 
for this Hamiltonian density has the general form 



where Oj are operators that depend on the chiral angle 9 of the soliton and 
the spin and radial coordinates of the mesons. In the matrix element (2.29) 
the meson fields M are replaced by the corresponding wave functions, which 
in the case of both mesons being in the ground state have the form (2.27). 
The general expression for the anharmonic correction (2.29) will then be 



Note that in the case when both mesons are in the same orbital symmetry 
they form a spin triplet state because they satisfy Bose statistics. 



3. The anharmonic correction for cascade hyperons 

In the case of two-meson states, when both mesons have the same flavour 
and are in the ground state (2.27) the expression for the matrix element of 
the anharmonic Hamiltonian derived above can be reduced to a very compact 
form, as in this case f ■ L = —2 and a 1 ■ a 1 = 1 as the mesons form a triplet 
state. The anharmonic correction that arises from the quadratic term in the 
Hamiltonian, which is obtained by converting C^f> (2.17) into a corresponding 
Hamiltonian density is 




(2.29) 



n (4) = M() : M.\r()- 2 .\i. 



(2.30) 




3irx 2 f 2 



1 



/ 



drr¥(r){ W 2 + i 9 (r)-^}, 



(3.1) 
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where the function g{r) is defined as 



g(r) = cosQcos 2 -. 



(3.2) 



The correction to the energy that arises from the quartic stabilizing term 
(2.20) takes the following form: 



Finally the chiral symmetry breaking term C% SB (2.24) yields the correc- 
tion 



Note that all the coefficients of the wave functions in these expressions 
are nonsingular at the origin. 

Numerically we find the anharmonic quartic correction to the energy of a 
two-kaon state to be -26 MeV (Table I). This represents only a 7% increase 
of the energy 390 MeV that is obtained in the harmonic approximation for 
the two meson state. The relative smallness of the anharmonic correction in- 
dicates that the harmonic approximation is sufficiently reliable to use in the 
study of hyperon structure and that thus the fairly satisfactory prediction 
for the spectra [4] and magnetic moments [18-20] of the strange hyperons 
have uncertainties of no more than 10% due to the anharmonic terms. Note 





-6rV 2 (r)]}- 



(3.3) 



A4 4 L = -3^2 / drr 2 k 2 (r){m 2 k 2 (r) 
-(l-^)[A: 2 (r)[-^ + V 2 + -A,(r)] 



(3.4) 
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that the 10% correction in the meson energy corresponds to only a 2% cor- 
rection to the predicted energies of the cascade particles, and therefore may 
be viewed as insignificant. The contributions AE^ 4 \ and AE^ SB are 

listed separately in Table I. Of the value -32 MeV for AE^ SB the mass term 
alone contributes -29 MeV, and it thus represents the dominant term. 

The anharmonic correction to the energy of a two .D-meson state is -233 
MeV (Table I). As the energy of the bound two D-meson state is 2452 MeV 
this represents a correction of only 10 percent. In the case of the 5-meson 
we find with x — 2.8 the bound state energy 7056 MeV. In this case the 
anharmonic correction to be -639 MeV and transfer less significant than in 
the case of the charmed cascade hyperon. The separate contributions to the 
total anharmonic correction are given in Table I. 

The fact that the relative magnitude of the anharmonic correction does 
not increase with increasing meson mass may be understood with the help of 
the following argument. At large distances the meson wave function behaves 
as 



where AMs a constant. In the case of the ground state (2.27) k(r) is finite 
and nonzero at the origin. The following approximate form for the meson 
wave function 



has the correct behaviour at the origin and exponential tail at large distances. 
The normalization constant N is determined by the normalization condition 



k(r) = N 




(3.5) 



r 



k{r) ~ Ne W '' 



(3.6) 



[1] 




For the wave function model (3.6) this implies that 



N ~ 



(m 2 - a; 2 ) 3 / 4 



(3.7) 
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when the meson mass and energy is large (uj » A(0)). For large energies the 
main anharmonic contribution is obtained from those terms in (3.1), (3.3) 
and (3.4) which are proportional to m 2 (the mass term) and uj 2 . If the soliton 
functions these terms in the integrands are replaced by constants they yield 
the estimate 

/V 4 i ? 

Here we have used the fact that uj ~ m for large meson mass values. By 
(3.7) it then follows that 

AE «, „ K-f^. (3.9, 

In ref. [14] it has been shown that the consistency of the soliton model 
requires that 

f~yfr, (3-10) 

and in ref. [21] that to a good approximation for large meson mass values 
the meson energy is given as 

co = ^=. (3.11) 

Combination of these results yield the following scaling law for the anhar- 
monic correction: 

AE ]=, (3.12) 



which shows that the anharmonic correction falls with meson mass and would 
become insignificant in the case of very large meson mass values. 



4. Discussion 



The results presented above for the anharmonic correction to the energy 
of a bound two-meson state show that it reduces the predicted energy in the 
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harmonic approximation for a two meson state by 10 %. This indicates that 
the harmonic approximation is reliable for the strange, charm and bottom 
hyperon generations, and that the satisfactory results it leads to for the static 
properties of the heavy flavour hyperons are fairly stable. As the topological 
model implies a large colour number approximation its inherent accuracy 
should in any case not be expected to be better than about 10%. It should 
also be noted that small adjustments of the values of the decay constants 
that appear in the Lagrangian density can lead to shifts in the predicted 
meson energies that are much larger than 10% [3]. 

The numerical results in this work apply only to the model with the alter- 
nate quartic stabilizing term (2.4). But the qualitative conclusion that the 
anharmonic corrections to the hyperon energies are small is not restricted 
to this particular version of the soliton model. It has been shown [14] that 
models with alternate stabilizing terms give qualitatively similar results in 
the harmonic approximation. The consistency argument for the smallness of 
the anharmonic correction given in section 3 is also general. 

The anharmonic terms do also in principle contribute to the energy of a 
single meson state in the form of self energy corrections. These corrections, 
which have not been considered here, require summation over all meson states 
and hence a method of regularizing the high energy terms in the sum. As the 
Lagrangian model (2.1) is constructed as an effective low energy interaction 
the evaluation of these self energy corrections lead outside the range of its 
validity and would require a systematic approach along the lines of chiral 
perturbation theory. 
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Table I 



Numerical values (in MeV) for the different anharmonic contributions to 
the energy of a 2-meson state with the mesons in the ground state. 





K 


D 


B 




3 


-6 


-21 


a4 4) 


3 


-44 


144 


a4 4 L 


-32 


-183 


-762 




-26 


-233 


-639 
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